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We compare the extraction of the ground-state decay constant from the two-point correlator in 
QCD and in potential models and show that the results obtained at each step of the extraction 
procedure follow a very similar pattern. We prove that allowing for a Borel-parameter-dependent 
' effective continuum threshold yields two essential improvements compared to employing a Borel- 

' parameter-independent quantity: (i) It reduces considerably the (unphysical) dependence of the 

, extracted bound-state mass and the decay constant on the Borel parameter, (ii) In a potential 

04 ' model, where the actual value of the decay constant is known from the Schrodinger equation, a 

^ ' Borel-parameter-dependent threshold leads to an improvement of the accuracy of the extraction 

, procedure. Our findings suggest that in QCD a Borel-parameter dependent threshold leads to a 

■ more reliable and accurate determination of bound-state characteristics by the method of sum rules. 

\0 ' PACS numbers: 11.55.Hx, 12.38.Lg, 03.65.Ge 
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INTRODUCTION 



^ In a series of recent publications we studied the extraction of the ground-state parameters from SVZ sum rules [l[ 

^ ■ (see also, e.g., 3)- We made use of a quantum-mechanical potential model since this is essentially the only case 
, where the standard procedures adopted in the method of sum rules may be tested: the estimates for the ground-state 
'~~'> parameters obtained by these procedures may be compared with the actual values of the ground-state parameters 

calculated from the Schrodinger equation, thus providing an unambiguous check of the reliability of the method. 
^ . The main results of our papers may be summarized as follows: (i) The standard approximation of a constant 
[ — ' effective continuum threshold does not allow one to probe the accuracy of the extracted hadron parameter [1, 0, 0, [1] ■ 
T-H (ii) Allowing for a Borel-parameter-dependent effective continuum threshold (we denote the Borel parameter r in 
QCD and T in the potential model) and fixing this quantity by using the information on the ground-state mass leads 
to a considerable improvement of the accuracy of the method [9] . 
PsJ The goal of this letter is to demonstrate that the results obtained at each step of the extraction procedure both in 

QCD and in potential models follow the same pattern. This similarity gives a strong argument that all our findings 
concerning the extraction of bound-state parameters from correlators obtained in potential model apply also to QCD. 
In particular, it points out a way of improving the results for bound-state parameters obtained from various correlators 
in QCD. 

The paper is organized as follows: In Section[2]we recall the QCD results from Ref. [l^] for the vacuum correlator of 
two pseudoscalar currents — the basic object for the extraction of fs within the framework of QCD sum rules. Section 
[3] provides the analogous results for a quantum-mechanical model for the case of a potential containing confining and 
Coulomb interactions. Section 2] compares the procedures of extracting the decay constant. Section O summarizes our 
conclusions. 

2. CORRELATOR AND SUM RULE IN QCD 

Let us consider the correlator 

n(p2) ^ ^ I d^xe^P^{0\T [Mx)jI{0)) |0) (2.1) 

of two pseudoscalar currents j^ix) — {mi, + mu)q{x)i^5b{x). The Borel-transformed operator product expansion 
(OPE) series for this correlator has the form 

n(r) = J e"'''^/9port(s,Ai)ds-|-npowcr(T,/i), (2.2) 

(niij-'r'mu)'^ 
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where the perturbative spectral density reads 

P,.rt{s, = m) + ^P^^H^, + f P^'^ {s,p) + ---, (2.3) 



/i being the renormalization scale. We make use of the results for pport reported in [T^l and do not reproduce the 
explicit expression for this quantity here. Following the argument of jlO| we work in terms of the running masses 
in the MS scheme. Therefore, in all expressions in this section the quark masses and to„, and as are the MS 
running quantities at the scale fi. Recall that the full Borel-transformed correlator (|2.2p does not depend on the 
renormalization scale /i; however, both the perturbative expansion truncated to a fixed order in as and the truncated 
power corrections depend on /i. We provide numerical estimates for /i — mi,; for this choice of the scale the known terms 
of the perturbative expansion exhibit a good hierarchy. We set mi,{mi,) — 4.2 GeV and for other QCD parameters 
make use of the central values reported in Table I of lO].^ The power corrections have been also considered in |10|]: 

npowcr(T, fi = mb) = {mb + m^J^e"™'^ (2.4) 



-mb{qq) 



2CFas A mlT\ 2 ^'^u ml ^ .2 



The parameter mg describes the contribution of the four-quark condensate. Notice that radiative corrections to the 
condensates increase rather fast with t. 

The correlator (|2.ip may be calculated in terms of hadron intermediate states: 

n(T) = ng(r) + contributions of excited states, ng(r) = flMge^^'^^'' , (2.5) 

where Jb is the decay constant of the i3-meson, defined by 

{mb + mu){0\mj5b\B) = /bM|. (2.6) 

For large values of r the contributions of the excited states decrease faster than the ground-state contribution and 
n(T) is dominated by the ground state. Unfortunately, the truncated OPE does not allow to evaluate the correlator at 
sufficiently large r, so the excited states give a sizeable contribution to II{t) for the considered values of r. 

According to the duality assumption, the contribution of the excited states is described by the perturbative contri- 
bution above some effective continuum threshold Scff . Then one obtains the following relation: 



ng(T) = ndual(T,Seff) (2.7) 



with 



Soft(T") 

ndual(T, Soff) = J e'''''ppeTt{s,p)ds + npowcr{T,p). (2.8) 

In the region near the physical continuum threshold at s = (Mb* + TO,r)^, the perturbative spectral density and the 
hadron spectral density are very different. Consequently, the effective continuum threshold as defined in (j2.7p turns 
out to be necessarily a function of the Borel parameter r. 

The necessity of the r-dependence of See may be understood by comparing the left-hand side (l.h.s.) and right-hand 
side (r.h.s.) of (|2.7p : the only way to obtain a single exponential on the l.h.s. for a given spectral density of the integral 
representation on the r.h.s. is to have a r-dependent Soff- 

The T-dependence of s^e may be also demonstrated explicitly: for any value of the ground-state parameters on the 
l.h.s. of (|2.7p one can obtain numerically Soff and see that it does depend on r. 

One should be aware of the fact that the r-dependence of Scff cannot and does not contradict to any principles of 
field theory: the dual correlator is a hand-made object; such an object does not emerge in field theory. Therefore, the 



It is well known that the numerical value of the correlator depends sizeably on the values of ni(,(mf,) and on the specific choice of the 
renormalization scale fi. However, a discussion of this dependence is far beyond the scope of this paper. A detailed analysis of fs in 
QCD is deferred to a separate publication. 
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properties of the dual correlator (e.g., its analytic properties) are very different from the properties of the field-theoretic 
correlators. 

Clearly, the standard assumption of a r-independent Seff is a possible assumption. We shall demonstrate, however, 
that relaxing this assumption leads to a visible improvement of the obtained results. 
We define the dual decay constant and the dual invariant mass by the relations 

fluAr) = M^4e*^^-nduai(T, Seff (r)), MiUr) - - 1: log Hduai (r, s,s{r)). (2.9) 

Notice that the deviation of the dual mass from the actual mass of the ground state gives an indication of the 
excited-state contributions picked up by the dual correlator. 



3. CORRELATOR AND SUM RULE IN POTENTIAL MODELS 



In parallel to QCD, let us consider a quantum-mechanical model with a potential containing a confining part, for 
which we take the HO form, and an attractive Coulomb interaction: 



H 



2m 



a 
r 



A quantum-mechanical analogue of the Borelized two-point function has the form ^ 

n(r) = (t^ = 0| exp(-i7T)|?^== 0). 



(3.1) 



(3.2) 



We construct the analogue of the OPE series for this correlator by retaining, similar to the QCD case, the perturbative 
contributions up to O(a^) (three loops of the non-relativistic field theory) and two power corrections, including 0{a) 
corrections to them. The resulting expression reads (the correlator for a pure Coulomb potential can be found in ll|):^ 

noPE(T) = npert(T) + npowcr(r), 



npcrt(T) 
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(3.3) 



Now, according to the standard procedures of the method of sum rules, the dual correlator is obtained as follows: 
we represent the perturbative contribution as a single spectral representation in the relative kinetic energy z of the 
interacting quarks and cut this representation at Zcfr: 



/ TO \ "^/"^ f I Z 

nduai(r, Zeff) = j / dzexp(-zT) 2W- 



'-+V2T 



3V^ 



n 



power 



(T). 



By construction, the dual correlator is related to the ground-state contribution by 

nduai(T,Zeff) = ng(T) = i?gexp(-£;gT), i?g = |V'g(r = 0)1^. 



(3.4) 



(3.5) 



As we have shown in our previous studies of potential models, the effective continuum threshold defined according to 
p.5[) is a function of the Borel time parameter T. 

For our numerical analysis, we adopt the following parameter values: a reduced quark mass of to = 0.175 GeV, which 
corresponds to a constituent quark mass of 0.350 GeV relevant for nonrelativistic computations; lo — 0.5 GeV, which 
leads to a realistic radius of the qq system, and a — 0.3. 

The energy and wave function of the ground state are found by solving numerically the Schrodinger equation with the 
help of the Mathematica code provided in ^]: the ground-state energy is Eg — 0.6473 GeV [for comparison, a pure HO 
model yields E^'^ — 0.75 GeV]; the ground-state wave function at the origin has the value V'(r = 0) 0.0783 GeV^/^ 
j^HO^j, — Q~) — (muj/'n:Y/^ = 0.068 GeV^^^]. Obviously, the effect of the Coulomb interaction is not small. 



^ Interestingly, at small values of T the system behaves like a free system, since the contribution of the confining potential as well as that 
of the radiative corrections vanish for small T. According to Q this is a signature of asymptotic freedom. So, a non-relativistic potential 
model 1 13.11 1 even with a constant a behaves like an asymptotically free theory. 
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4. EXTRACTION OF THE DECAY CONSTANT 

Whether a r-independent or some r-dependent effective continuum threshold is considered, the crucial problem is the 
choice of the criterion for fixing this quantity. We shall proceed as follows: 

A. The Borel window 

First, we must fix the working r window where, on the one hand, the OPE gives an accurate description of the exact 
correlator (i.e., the higher-order radiative and power corrections are small) and, on the other hand, the ground state 
gives a sizeable contribution to the correlator. In QCD, we set the window as follows: 

0.05 GeV"2 < r < 0.18 GeV~^ (4.1) 

In the region r < 0.18 GeV~^ the as and a1 terms to n(r) contribute less than 10% and 3% of the leading term, 
respectively. Power corrections give about 20% of the leading term. We point out that the radiative corrections to 
the condensates increase rather fast with r, so it is preferable to stay at relatively low values of r. Therefore our 
window is located at the lower values of r compared to the window adopted in fl^ ; in the region (|4.1|) the accuracy 
of the truncated OPE is higher. 

It is known that the experimental value of the B meson decay constant is fs ~ 200 MeV. Adopting this value, we 
may calculate the relative contribution of the ground state to the correlator. In our window it does not exceed 50%. 

In the potential model, we choose the window as 

0.2 GeV"^ < T < 0.8 GeV"^ (4.2) 

For these values of T the omitted unknown higher-order power corrections are negligible, so the correlator is known 
with good accuracy. The relative contribution of the ground state to the correlator amounts to 10% at T = 0.2 GeV^^ 
and to 50% at T = 0.8 GeV'^ 

We shall see that for a relative ground-state contribution of this size our procedure allows one to extract the decay 
constant with a reasonable accuracy. 

B. Fixing the effective continuum threshold 

Widely used is the so-called stability criterion: one looks for that constant value of Scs for which the extracted decay 
constant is most stable in the window. This Borel stability is an implementation of a self-evident statement that the 
physical observable cannot depend on r, an auxiliary parameter of the method. The problem is, however, that the 
independence of a hadron decay constant of r, being a necessary condition, is not sufficient to guarantee the extraction 
of the right value. We have given several examples for potential models which nicely demonstrate that assuming a 
T-independent effective continuum threshold and fixing its value by requiring maximal stability in the Borel window 
may lead to the extraction of a very inaccurate value. 

In this paper, we consider a different algorithm for the extraction of /b, which makes use of the knowledge of the 
ground-state mass [l3|- In parallel to QCD we present also the results for a quantum- mechanical model p.ip . This 
is done in order to demonstrate the way our algorithm works in the case where the exact value of the decay constant 
is known. A comparison makes clear that, with respect to the extraction procedure, there are no essential differences 
between QCD and quantum mechanics. 

The algorithm developed in our previous works and established to work well for different correlators in the potential 
model is very simple: we consider a set of r-dependent Ansatze for the effective continuum threshold (for the case of 
the potential model one just replaces r ^ T and M ^ E): 

n 

4ff^W-E4"^-^- (4-3) 
i=o 

Obviously, the standard r-independent effective continuum threshold is also taken account by (|4.3p . Now, we fix the 
parameters on the r.h.s. of (|4.3p as follows: we calculate the dual mass squared according to (|2.9p for the r-dependent 
Seff of Eq. (|4.3p . We then evaluate M1^^^{t) at several values of r = ri (i = 1, iV, where N can be taken arbitrarily 
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Fig. 1: Left column: potential model (13. If) : right column: QCD. First line: relative contribution of the ground state to the 
correlator; second line: fitted dual mass; third line: corresponding dual decay constant. The dashed line in Fig. (e) corresponds 
to the true value of the decay constant obtained by solving the Schrodinger equation. The index n is the power of the polynomial 
Ansatz for the Borel-parameter-dependent effective continuum threshold. 



large) chosen uniformly in the window. Finally, we minimize the squared difference between AfJ^^j and the known 
value Mg-. 

x'^^E[^^dual(-.)-M|]^ (4.4) 

1=1 

This gives us the parameters of the effective continuum thresholds. As soon as the latter are fixed, it is straightforward 
to obtain the decay constant. 

Figured] shows the results for QCD and for our potential model; in the latter the actual value of the decay constant 
has been found from the Schrodinger equation, so that we may control each step of the extraction procedure. 

First, let us notice that the Borel-parameter-dependent effective thresholds corresponding to n = 1 and n = 2 lead 
to a visible improvement of the stability of the dual mass [Figs. [TJc) and (d)] compared to the constant threshold. 
This means that the dual correlator for n > is less contaminated by the excited states; according to the philosophy 
of QCD sum rules the better stability of Afduai with ?i > is an important achievement for the trustability of the 
results. 

According to Fig. [U in the potential model the true value of the decay constant lies in the band provided by the 
linear (n — 1) and the quadratic (n = 2) fits. We have checked that this result holds in a broad range of the parameters 
of the potential model. The similarities of each step of the extraction procedure in QCD and in the potential model are 
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evident.^ Therefore, it is tempting to expect that also in QCD the decay constant hes in the range provided by the 
hnear and the quadratic fits. Anyway, the difference of the results obtained for n — I and n = 2 constitutes a realistic 
estimate of the intrinsic uncertainty of the extracted decay constant. If one considers only the standard constant 
Ansatz (n = 0) for the effective continuum threshold, the accuracy of the extracted decay constant cannot be probed. 

5. DISCUSSION AND CONCLUSIONS 

We have presented a detailed analysis of the extraction of the decay constant from the two-point function in QCD 
and in a potential model. Our results may be summarized as follows: 

(i) The comparison presented in this work makes obvious that, with respect to the extraction of the ground-state 
parameters, there are no essential differences between QCD and quantum mechanics: as soon as the parameters 
of the Lagrangian are fixed, and the truncated OPE is calculated with a reasonable accuracy (taking into account 
also the relevant choice of the renormalization scale in QCD), the extraction procedures are very similar. 

At first glance, this similarity might look surprising since we know that the structure of bound states in potential 
models and in QCD are rather different. However, the method of dispersive sum rules does not make use of 
(and does not provide) information about the details of the ground-state structure. What really matters for 
extracting the ground-state parameters in this method is the structure of the OPE for a given correlator. Since 
the structure of the OPE in QCD and of its analogue in potential models is rather similar, it should not be 
surprising at all that the extraction procedures in potential models and in QCD are similar, too. 

In view of the above similarity, our previous results for the extraction of the ground-state parameters (including 
also the form factors [9]) obtained in potential models have direct implications for the corresponding analyses in 
QCD and should be taken quite seriously. 

(ii) Allowing for r-dependent Ansatze for the effective continuum threshold leads to two essential improvements: 

(a) The stability of both the dual mass and the dual decay constant in the window is considerably improved if 
one proceeds from the standard r-independent to the r-dependent Ansatz for the effective continuum threshold. 

(b) In the potential model, where the exact decay constant has been calculated from the Schrodinger equation, 
allowing for a r-dependent effective continuum threshold and fixing its parameters according to (j4.4|) — i.e., 
by minimizing the deviation of the dual mass from the known ground-state mass in the window — leads to 
the extraction of a more accurate value. As follows from our analysis, a realistic band of values of the decay 
constant is provided by the numerical results obtained with the linear and quadratic Ansatze for the effective 
continuum threshold. The intrinsic uncertainty (i.e., the one related to the extraction procedure) of the decay 
constant found in this way is expected to be at the level of a few percent. 

Although not rigorous in the mathematical sense, this estimate for the systematic uncertainty may be considered 
as a realistic educated guess supported by findings in models where the true value of the decay constant is known. 
Moreover, we have doubts that a more rigorous estimate of the intrinsic error of the method of sum rules may 
be obtained in principle. 
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^ One may observe that outside the window the behavior of the dual mass in QCD and in the potential model is not exactly the same. 
This is related to the fact that the quark condensate in QCD is negative, whereas the corresponding power correction in potential models 
(for any confining potential) has a positive sign. 
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